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Abstract 

An effective mass scale for supersymmetric particles is defined and techniques 

for its measurement at the LHC discussed. Monte Carlo results show that, for jets + 
^miss events, a variable constructed from the scalar sum of the transverse momenta of all 
reconstructed jets together with E^^'^'^ provides in most cases the most accurate model 
independent measurement of MJ'^y (intrinsic precision ^ 2.1 % for mSUGRA models). 
The overall precision with which M^^^y could be measured after given periods of LHC 
running and for given classes of SUSY models is calculated. The technique is extended to 
measurements of the total SUSY particle production cross section cTsusy 
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1 Introduction 



One of the principal motivations for construction of the Large Hadron Colhder is the search 
for low energy supersymmetry (SUSY) Q. In a large class of models the interactions of SUSY 
particles conserve R-parity, causing the Lightest Supersymmetric Particle (LSP) to be neutral 
and stable. R-Parity conserving SUSY events at hadron colliders are predicted to consist 
of cascade decays of heavy, strongly interacting SUSY particles into lighter Standard Model 
(SM) particles and two LSPs. This results in the classic discovery signature of an excess of 
events containing jets, leptons and large quantities of event missing transverse energy E!^^^^ 

i- 

Should R-Parity conserving SUSY particles be discovered at the LHC the next task would 
be to measure their properties. Importantly, such measurements must be independent of the 
SUSY model and its parameters, which are a priori unknown. This process is complicated by 
lack of knowledge of the momenta of the two escaping LSPs in each event, preventing direct 
reconstruction of SUSY particle masses. Consequently other techniques are required which can 
measure masses or combinations of masses indirectly. With sufficient integrated luminosity it 
should be possible to look for edges in the invariant mass spectra of various combinations of 
jets and leptons in SUSY events |2|, but initially the most effective technique is likely to be 
the use of distributions of event transverse momentum px and missing transverse energy Q . 
In this letter, the latter technique will be investigated in detail and extended to SUSY models 
beyond the minimal Supergravity (mSUGRA) models considered previously ||3|. 



2 Measurement Technique 

Consider a heavy SUSY particle (mass mi) produced in a hadron-hadron collision. Assume 
further that this particle is boosted along the beam-axis by the longitudinal momentum im- 
balance of the event. If this particle undergoes a cascade decay to a lighter SUSY particle 
(mass 7712) and a Standard Model particle (assumed massless), then the transverse momentum 
Pt of the SM particle in the lab frame is related to nii and m2- 

1 / 

PT cc-lmi ^ . (1) 

2 V mi J 

Variables based on the px of SM particles in SUSY events are therefore sensitive to SUSY 
particle masses, modulo smearing effects arising from the true rj distribution of those particles. 

Events in the jets -|- E!^^^^ + leptons channel were used for this study. The lepton 
veto requirement was imposed to reduce possible systematics in the measurement arising from 
SM neutrino production. Defining PT{i) as the transverse momentum of jet i (arranged in 
descending order oi pt), the following four measurement variables Mest were studied: 

(1) Mest = \PT(1)\ + \PT{2)\ + IPT{3)I + IPT{4)I + Ef^^', 

(2) Mest = |PT(1)I + \PT{2) \ + \PT{3) \ + IPT(4)|, 

(3) Mest = E^|PTwl+i?^^ 

(4) Mest = E^|PT{^)|■ 
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The first variable is identical to the "effective mass" variable (-Meg) defined in Ref. Q. The 
scalar sum of the transverse momenta of only the four hardest jets was used due to the 
predominantly four-jet nature of many SUSY events [Q], while the addition of the event E^^^^ 
accounts for the pT carried away by the LSPs. 

Combining all jets in each detector hemisphere j into one pseudo-particle of transverse 
momentum PT{j) and invariant mass mj, a fifth variable was also defined: 

(5) Mest = i E?=i ^rnj + 4p^^^.^ + 2^ {ml + 2pl^^^).{m] + 2pl^^^), 

where rrix = 100 GeV/c^. This variable approximates the mean of the reconstructed masses 
of the two initial SUSY particles, assuming that each moves close to the beam axis and decays 
into an LSP with mass of order rrix. 



3 Definition of Mass Scale 

Due to the large number of different SUSY particles which can be produced in any given event 
masses measured with these variables will not correspond to those of any one particular SUSY 
state. Nevertheless in many models the strongly interacting SUSY particles are considerably 
heavier than states further down the decay chain. Hence it is the decays of these particles 
which will contribute most to the sum of \pt\- For this reason we shall choose to define a 
SUSY "mass scale" Mgusy as the weighted mean of the masses of the initial SUSY particles, 
with the weighting provided by the production cross section of each state: 



M,„3y = (2) 



This definition differs from that used in Ref. Q but in the limit where squarks or gluinos of 
a single mass dominate the production cross section it gives the same result. 

The above approximation breaks down for models where the lighter SUSY particles are 
of similar mass to the strongly interacting states. We shall attempt to compensate for this 
when using variables (1) - (4) by defining an effective SUSY mass scale M^^^y in analogy with 
Eqn. (§: 

= {"■•"' - • 

where is the mass of the LSP. For variable (5) the equivalent expression is somewhat 
different: 



^sfsy = V^slsy-M2. (4) 

With these definitions of the effective SUSY mass scale comparison of experimental results with 
the predictions of a given SUSY model requires knowledge of the model dependent particle 
mass spectrum and production cross sections. 
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4 Simulation and Event Selection 

Events were generated using PYTHIA 6.115 g (SM background, mSUGRA and MSSM sig- 
nal) and ISAJET 7.44 (GMSB signal). Hadronized events were passed through a simple 
simulation of a generic LHC detector. The calorimeter was assumed to have granularity 
Arj X A(p = 0.1 X 0.1 over the range |r/| < 5, and energy resolutions 10%/^/E © 1% (ECAL), 
50%/VE © 3% (HCAL) and 100%/V^ © 7% (FCAL; |r/| > 3). Jets were found with the 
GET JET g fixed cone algorithm with cone radius AR = 0.5 and E!j^^ = 50 GeV. 
Events in the jets + E^^^^ channel were selected with the following criteria: 

• > 4 jets with pr > 50 GeV 

• > 2 jets with pr > 100 GeV 

• Sf?'^^ > max(100 GeY,0.25j2iPT{i)) 

• Transverse Sphericity St > 0.2 

• A'^{PT(.„PT(.,) < 170° 

• ^'^(PT(.,+PT,.),Er=) < 90° 

• No muons or isolated electrons with px > 10 GeV in |r/| < 2.5. 

Standard Model background events were generated for the following processes: tt (5 x 10^ 
events), W + jet (5 x 10^ events), Z + jet (5 x 10^ events) and QCD 2^2 processes Q 
(2.5 X 10*^ events). The distribut ions of the M^st variables for these events were then compared 
with those for SUSY signal events generated from the mSUGRA (minimal Supergravity ||^]), 
MSSM (Minimal Super symmetric Standard ModelQH) and GMSB (Gauge Mediated SUSY 
Breaking Q) models. In each case 100 points were randomly chosen from within the parameter 
space of the model, and 1 x 10^ events (a factor 10 greater than in Ref. Q) generated for each 
point. 

For mSUGRA models the region of parameter space sampled was identical to that used 
in Ref. g: 100 GeV < mo < 500 GeV, 100 GeV < mi/2 < 500 GeV, -500 GeV < Aq < 500 
GeV, 1.8 < tan(/3) < 12.0 and sign(^) = ib 1. For MSSM models the choice of points was 
complicated by the requirement that models be physically realistic (positive particle masses 
etc.) and be suitable for this study (neutralino LSP). The masses of the strongly interacting 
SUSY particles and sleptons were constrained to lie in the range from 250 GeV/c^ to 2000 
GeV/c^ while the mass parameters of the partners of the electroweak gauge bosons were 
constrained to lie in the range from 50 GeV/c^ to the mass of the lightest strongly interacting 
SUSY particle or slepton. tan(/?) was constrained to lie in the range 1.4 < tan(/3) < 100.0. 
Finally, given that the lightest SUSY particles in models with low Mg|^y have a high probability 
of discovery before LHC comes on line, only those models with M^^^y > 250 GeV/c^ were used. 

In mSUGRA and MSSM models the LSP is generally the lightest neutralino x?) but in 
GMSB models this role is taken by the gravitino. If the Next-to-Lightest Supersymmetric 
Particle (NLSP) is neutral and sufficiently long-lived to escape from the detector then the 

constrained version of the MSSM with 15 free parameters (no additional D-terms, 3rd generation trilinear 
couplings derived from masses) implemented in SPYTHIA 
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phenomenology is similar to that for mSUGRA or MSSM models Q. If the NLSP is short- 
lived however then it decays to a gravitino and the phenomenology is different. To test whether 
the Mest variables defined above are also sensitive to the effective mass scales of these latter 
models points were chosen from within the range of GMSB parameter space defined by: 10 
TeV < Am < 100 TeV, 100 TeV < M„ < 1000 TeV/c^, 1 < A^g < 5, 1.8 < tan(/3) < 12.0 
and sign(/x) = it 1. The value of Cgrav, the ratio of the gravitino mass to that expected for 
only one SUSY breaking scale |6|, was set to unity in all cases to ensure rapid decays to the 
gravitino LSP. Again only those models with Mg^y > 250 GeV/c^ were used. 

In the mSUGRA and MSSM models a statistically significant excess of signal events (S) 
above background (B) S + B - ^/B > 5.0 lol) was found for the majority of points after 
the delivery of only 10 fb^^ of integrated luminosity (1 year of low luminosity operation). In 
GMSB models the data indicate that greater event statistics corresponding to at least 100 fb~^ 
(1 year of high luminosity running) would be required for discovery with these particular cuts. 
It should be noted however that the use of this channel and these cuts has been optimised for 
mSUGRA points. In GMSB models with prompt decays to gravitino LSPs photon production 
(bino NLSP) or lepton production (slepton NLSP) is common [11] and consequently many 
events were rejected by the lepton veto and jet multiplicity requirements. In dedicated GMSB 
studies these requirements should therefore be loosened in order to increase signal acceptance. 
In this case measurement variables taking account of lepton and/or photon px should also be 
used to reduce systematic measurement errors |11]. 



5 Mass Scale Measurement 

The Mest distributions of SUSY signal events in the models considered here are roughly gaus- 
sian in shape (see Figs. (1) - (5) of Ref. Q), in sharp contrast to the SM background which 
falls rapidly with Mest. Fitting gaussian curves to the signal distributions then provides es- 
timates of their means which can be compared with the effective mass scales M^i^y of the 
corresponding SUSY models. The degree of correlation between the two variables gives a 
measure of the intrinsic (systematic) precision provided by the Most variable when measuring 

A typical scatter plot of M^l^^^ against Mgst (here variable (3)) for mSUGRA models is 
shown in Fig. ||(a). The correlation between the variables is clearly very good. To quantify the 
degree of correlation a linear regression was performed on the data and the points projected 
onto an axis perpendicular to the fitted trendline. The distribution of the data along this line 
for mSUGRA models is shown in Fig. |2|(a). 

A scatter plot of M^^^^ against Afcst variable (3) for MSSM models is shown in Fig. |l|(b). 
In Figs. |2|(b) and (c) are plotted correlation histograms derived from this figure using the pro- 
jection axis defined by the mSUGRA data. Fig. |2|(b) shows the histogram obtained assuming 
-^sSy = Msusy, while Fig. ^(c) shows the equivalent histogram using the definition of Mg°^y in 
Eqn. (^). The smaller scatter in this latter case indicates an improved measurement precision. 
An improvement is also obtained for mSUGRA models. Here however it is smaller since is 
usually much less than the masses of the strongly interacting SUSY particles ||l^. In GMSB 
models (Fig. |(c) and Fig. |(d)) M^ is negligible and so the two definitions of Mg'^^y are always 
identical. 

The correlation histograms were next fitted with gaussian functions. The fitted values 
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widths a were used to calculate the intrinsic measurement precision by subtracting in quadra- 
ture the width expected from finite event statistics (assumed to be given by the rms of the 
errors on the fitted means of the Mest distributions). Due to uncertainties in the expected 
statistical scatter this correction could give unduly optimistic estimates of the measurement 
precision if it were large. In all cases however the corrections were found to be small (< 33%) 
and the effects of such uncertainties neglected. 

Intrinsic measurement precisions for M|^y calculated using the above technique for mSUGRA, 
MSSM and GMSB models and the five Mgst variables listed in Sec. 2 are presented in Table |l[ 
Variable (3) (Mest = SilPT(i)l + -E™'****) provides the greatest precision for mSUGRA and 
MSSM models, with the higher precision being for mSUGRA models (2.1 %). The poorer 
precision for MSSM models (12.8 %) is due to their greater number of free parameters and 
hence the smaller correlation between the particle masses. In the GMSB models variable (4) 
(Mest = X]j |PT{j)l) provides the greatest measurement precision (6.1 %), however variable (3) 
is reasonably accurate (9.0 %). This indicates that effective mass scale measurements are also 
effective for models with gravitino LSPs. 

It should be noted that for any given Mest variable the fitted means of the projected 
histograms (Table |) for mSUGRA, MSSM and GMSB models are consistent to within the 
fitted widths. For this reason it can be said that the variables provide model independent 
measurements of the effective SUSY mass scale M^^^y. The expected measurement precision 
is SUSY model dependent (due to larger widths for MSSM histograms than for mSUGRA and 
GMSB histograms) but this is less troublesome when comparing measurements with theory 
because in this case a particular SUSY model must be assumed. 

With the intrinsic precision from Table ^ for M^st variable (3) it is possible to estimate 
the overall (systematic + statistical) precision for measuring effective SUSY mass scales in 
mSUGRA, MSSM and GMSB models as a function of the mass scale and integrated luminosity. 
Distributions of Mgst for signal + background events were first constructed assuming integrated 
luminosities of 10 fb~^ (1 year low luminosity), 100 fb^^ (1 year high lumi.) and 1000 fb^^ (10 
years high lumi.). The mean background distribution was then subtracted from each with an 
assumed 50% systematic error 0. The resulting distributions were again fitted with gaussian 
functions and the errors on the fitted means added in quadrature to the intrinsic precision 
calculated above and an estimated 1% systematic error from uncertainties in the measurement 
of the jet energy scale. 

The results are plotted in Fig. ^ for mSUGRA, MSSM and GMSB models. In the last case 
results for 100 fb~^ and 1000 fb~^ only are presented due to the poor statistical significance 
of GMSB models in the jets + E!^^'^'^ channel at low integrated luminosity. Precisions < 
15 % (40 %) should be achievable in mSUGRA (MSSM) models after only one year of low 
luminosity running, improving to < 7 % (20 %) after one year of high luminosity running. 
Due to the poor statistics obtainable from GMSB models, particularly for high Mg|^y values, 
measurement precisions < 50 % are likely to be obtainable with these cuts only after the 
delivery of 1000 fb^^ of integrated luminosity and for M^l^ < 1000 GeV/c^. 



^It is unlikely that the distribution of background events (especially QCD) will be known with any certainty 
from theory. Instead the distribution of low events will likely be measured and the results extrapolated 

into the high iJ™'^'^ region [ ^ , 50% is a conservative estimate of the systematic uncertainty associated 
with this extrapolation. 



5 



6 Cross Section Measurement 



The above technique can also be used to measure the total SUSY production cross section 
(Tsusy, although in this case it is the fitted normalisation of the signal Mest distribution which is 
of interest. The correlation between this normalisation and Csusy for M^st variable (3) is shown 
in Fig. ^. The correlation is reasonably good with the data best fitted by a power-law. For 
these measurements the errors are non-gaussian (due to the power-law relation between the 
normalisation and Csusy) and in reality it is the logarithm of the measured cross sections which 
is approximately gaussian distributed. For this reason the intrinsic measurement precisions for 
In ((Tsusy) are listed in Table |. The overah (non-gaussian) precisions for CTgusy (i.e. dcTg^gy/cXgi^gy 
= d(ln((Tsugy))) are plotted in Fig. |5| in the same format as Fig. |^. For mSUGRA models the 
overall figugy measurement precision obtainable for 1000 fb~^ is < 15 %, while for MSSM 
models it is < 50 %. 

Measurements of cigugy carried out in this way are inherently sensitive to the type of 
SUSY model, in contrast to the measurements of Mjf^y. This is because in some models (e.g. 
GMSB) the SUSY particle decay characteristics can be such that the probability for signal 
events to pass the selection cuts is reduced significantly relative to that for mSUGRA models. 
The analysis presented here is intended to be model independent and so projects data onto a 
single axis perpendicular to the trendline of the mSUGRA models. Consequently in the GMSB 
case, where the trend is very different from that for mSUGRA, the presented measurement 
precision is poor (> 300 %). If it were known that GMSB models were correct then an axis 
perpendicular to the GMSB trendline could be used to obtain much greater measurement 
precision (ln(o"gugy) precision < 2.5 %). This highlights the fact that in reality measurements 
of fjgugy, unlike measurements of M^^^y, are dependent on the assumed SUSY model. 

7 Conclusions 

Model independent techniques for measuring the effective mass scale of SUSY particles at the 
LHC have been investigated. Overall measurement precisions better than 15 % (40 %) should 
be possible for mSUGRA (MSSM) models after only one year of running at low luminosity. 
Measurements should also be possible for models with rapid decays to gravitino LSPs, although 
with the requirement of either significantly increased statistics or measurement variables using 
photon or lepton pj-- The total SUSY production cross section should be measureable in a 
similar way ultimately to ~ 15 % (50 %) in mSUGRA (MSSM) models, although not in a 
completely model independent manner. 
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Tables 



Table 1: Estimates of the intrinsic M^^y measurement precision for mSUGRA, MSSM and 
GMSB models for the five M^st variables discussed in the text. The third and fourth columns 
show the fitted mean and width of the projected M^^^y - Mgst correlation histogram for each 
model and variable, and the fifth column their ratio. The sixth column shows the expected 
fractional width estimated from the rms error on the fitted means of the signal distributions. 
The seventh column contains the intrinsic measurement precision estimated by subtracting in 
quadrature column six from column five. 

Table 2: Estimates of the intrinsic In(cJsusy) measurement precision for mSUGRA, MSSM 
and GMSB models for M^st variable (3) (defined in the text). The third and fourth columns 
show the fitted mean and width of the projected ln(c7susy) - normalisation correlation his- 
togram for each model and variable, and the fifth column their ratio. The sixth column shows 
the expected fractional width estimated from the rms error on the fitted normalisations of 
the signal distributions. The seventh column contains the intrinsic measurement precision 
estimated by subtracting in quadrature column six from column five. 
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Figures 



Figure 1: The effective SUSY mass scale M|^y plotted against Mest for variable (3) (defined 
in the text) for 100 random mSUGRA (Fig. 1(a)), MSSM (Fig. 1(b)) and GMSB (Fig. 1(c)) 
models. Note the differing scale in Fig. 1(c) due to the larger spread in Mgst values generated 
for GMSB models. In Fig. 1(c) those GMSB models where the gaussian fit to the signal Mest 
distribution failed due to insufficient acceptance are omitted. 

Figure 2: Projections of the points in Fig. 1 onto an axis transverse to the fitted trend- 
line of mSUGRA data (Fig. 1(a)) for Mest variable (3). Fig. 2(a) shows the distribution for 
mSUGRA points, Fig. 2(b) the distribution for MSSM points with M^^^^ = M^nsj, Fig. 2(c) 
the distribution for MSSM points with M|^y given by Eqn. (^) and Fig. 2(d) the distribution 
for GMSB points. Bin widths are equal in Fig. 2(b) and Fig. 2(c) to aid comparison. Bin 
widths differ between other plots. 

Figure 3: Overall precision for measurement of M^^^^ after delivery of integrated luminosities 
of 10 fb~^ (stars), 100 fb~^ (open circles) and 1000 fb~^ (filled circles) for Mgst variable (3). 
Precisions for mSUGRA points are plotted in Fig. 3(a), MSSM points in Fig. 3(b) and GMSB 
points in Fig. 3(c). No data are shown for GMSB points for 10 fb~^ integrated luminosity due 
to the poor statistical significance of signal events in this scenario. Note the differing scale in 
Fig. 3(c) due to the larger spread in Mest values generated for GMSB models. 

Figure 4: The total SUSY particle production cross section (Tsusy plotted against the fitted 
normalisation of the signal distribution for variable (3) (defined in the text) for 100 random 
mSUGRA (Fig. 4(a)), MSSM (Fig. 4(b)) and GMSB (Fig. 4(c)) models. In Fig. 4(c) those 
GMSB models where the gaussian fit to the signal Mest distribution failed due to insufficient 
acceptance are omitted. 

Figure 5: Overall (non-gaussian) precision for measurement of cJsusy after delivery of integrated 
luminosities of 10 fb^^ (stars), 100 fb~^ (open circles) and 1000 fb~^ (filled circles) for Mest 
variable (3). Precisions for mSUGRA points are plotted in Fig. 5(a), MSSM points in Fig. 5(b) 
and GMSB points in Fig. 5(c). No data are shown for GMSB points for 10 fb~^ integrated 
luminosity due to the poor statistical significance of signal events in this scenario. Note the 
differing scale in Fig. 5(c) due to the larger spread in Mest values generated for GMSB models. 



9 



Model 


Variable 


X 


a 


a/x 


rms error (x) 


Precision (%) 


mSUGRA 


1 


1.585 


0.049 


0.031 


0.011 


2.9 




2 


0.991 


0.039 


0.039 


0.010 


3.8 




3 


1.700 


0.043 


0.026 


0.015 


2.1 




4 


1.089 


0.030 


0.028 


0.011 


2.5 




5 


1.168 


0.029 


0.025 


0.013 


2.1 


MSSM 


1 


1.657 


0.386 


0.233 


0.031 


23.1 




2 


0.998 


0.214 


0.215 


0.042 


21.1 




3 


1.722 


0.227 


0.132 


0.031 


12.8 




4 


1.092 


0.143 


0.131 


0.029 


12.8 




5 


1.156 


0.176 


0.152 


0.034 


14.8 


GMSB 


1 


1.660 


0.149 


0.090 


0.037 


8.1 




2 


1.095 


0.085 


0.077 


0.040 


6.6 




3 


1.832 


0.176 


0.096 


0.034 


9.0 




4 


1.235 


0.091 


0.074 


0.041 


6.1 




5 


1.273 


0.109 


0.086 


0.034 


7.9 



Table 1: 
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Model 


Variable 


X 


a 


a/x 


rms error (x) 


Precision (%) 


mSUGRA 


3 


0.855 


0.008 


0.009 


0.003 


0.8 


MSSM 


3 


0.848 


0.023 


0.027 


0.004 


2.7 


GMSB 


3 


0.742 


0.141 


0.190 


0.006 


19.0 



Table 2: 
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